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r9 : Abstract 



We consider a quasi one-dimensional chain of N chaotic scattering elements 

with periodic boundary conditions. The classical dynamics of this system is 

dominated by diffusion. The quantum theory, on the other hand, depends 

^^ I crucially on whether the chain is disordered or invariant under lattice transla- 

^ I tions. In the disordered case, the spectrum is dominated by Anderson local- 

^^ ■ ization whereas in the periodic case, the spectrum is arranged in bands. We 

r^ ■ investigate the special features in the spectral statistics for a periodic chain. 

For finite A'^, we define spectral form factors involving correlations both for 

identical and non-identical Bloch numbers. The short-time regime is treated 

'^ ' within the semiclassical approximation, where the spectral form factor can be 

Q ■ expressed in terms of a coarse-grained classical propagator which obeys a dif- 

2 ' fusion equation with periodic boundary conditions. In the long-time regime, 

Q . the form factor decays algebraically towards an asymptotic constant. In the 

limit iV ^ oo, we derive a universal scaling function for the form factor. The 

theory is supported by numerical results for quasi one-dimensional periodic 

;h ■ chains of coupled Sinai billiards. 
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I. INTRODUCTION 

The goal of quantum chaology is to identify universal features which prevail in systems 
whose classical dynamics is chaotic, and to relate the quantum observables to their classical 
counterparts. Spectral statistics is one of the topics where this approach has met with 
success, and the connection between the spectral form factor (the Fourier transform of 
the two-point correlation function) and the classical probability to stay, is a very fruitful 
junction of mesoscopic physics and quantum chaos |l]-§[. In the present article we develop 
the corresponding theory for extended, quasi one-dimensional systems. 

We consider composite systems consisting of A^ chaotic unit cells of length a with pe- 
riodic boundary conditions. The example we have in mind is a chain of interconnected 
chaotic billiards, as shown in Fig. |I]. In practice one can realize such systems by construct- 
ing a mesoscopic ring with a circumference L = Na which is smaller than the dephasing 
mean free path, but much larger than the elastic mean free path. The particle moves as a 
random walker between the unit cells, which results on average in a diffusive evolution. It 
is important to note from the outset that the diffusive nature of the classical evolution is 
completely indifferent to the translational symmetry of the ring — whether it is periodic or 
disordered. In contrast, quantum mechanics depends crucially on this feature: disordered 
rings are affected by Anderson localization and in the limit N ^ oo the spectrum remains 
discrete with a Poissonian spectral statistics 0. Chains which are invariant under discrete 
lattice translations develop continuous band spectra in the same limit. Can one account for 
these two different quantum situations within a theory that relies on the classical diffusive 
evolution which in itself is indifferent to the symmetry? The present work provides an af- 
firmative answer to this question. We study in particular the spectral form factor for the 
periodic case, and show that in the limit A^ -^ oo, it approaches a limit distribution which 
expresses the effect of the levels clustering into bands. 

Let us consider the two scenarios in some more detail: In the disordered case, once the 
overall length L = Na is sufficiently larger than the Anderson localization length ^, one can 
consider the system qualitatively as a union of L/^ uncoupled systems. Within this picture, 
the spectrum is a superposition of uncorrelated spectra, resulting in a spectral form factor 
Kj\f{T) ~ Ki{tL/C,), where the suffix 1 denotes the spectral form factor for the unit cell, as 
opposed to the form factor Kn{t) for the composite system. 

In the periodic case, the symmetry with respect to discrete lattice translations implies 
the existence of discretized energy bands. According to Bloch's theorem, the energies are 
labeled by a quasi-momentum (Bloch number) q which takes discrete values for a finite chain. 
Energy levels with different Bloch numbers are correlated over some range in q. Hence, on 
the time scale corresponding to an energy scale of the order of the inter-band spacing, the 
spectrum looks as if it were degenerate with a multiplicity ~ A^. This corresponds to a peak 
~ A^ in the form factor Km{t) at r ~ 1/A^. For larger times, finer structures in the spectrum 
are resolved, until for very large times (r ~ 1), the discrete nature of the spectrum is the 
dominant feature. On this time scale Kj\f{T) ~ 1 (in the absence of symmetries). As will be 
shown in this paper, the interpolation between the two time scales is given by a power-law 
decay K]^{t) ~ 1/r. 

Since we are dealing here with the periodic case, we make use of the extra parameter 
which quantum mechanics provides, and consider spectral statistics which test the corre- 



lations between spectra of different Bloch numbers. We develop a semiclassical theory for 
these spectral measures and show that they are intimately connected with the distribution 
of winding numbers in the symmetry-reduced classical system, which in the present case 
is the unit cell with periodic boundary conditions. This analysis makes use of the general 
formalism proposed by Robbins |0, and can be considered as an extension of the previous 
work on band spectra in chaotic systems |§]. 

It is important to emphasize that the considerations given above are rather general and 
do not depend on the details of the actual system. To check the universality of our arguments 
we have performed a parallel study where the system of interest is not a billiard chain, but 
an appropriately defined version of the quantum kicked rotor [Q . This work will be reported 
elsewhere |jlO|. The results obtained for the two models agree in all details, giving credence 
to our claim that the spectral statistics for periodic diffusive chains are universal. 

The paper is organized as follows. In section |I| we derive a secular equation for the 
energy spectrum of the billiard chain. It is based on the scattering approach to quantization 
|Tl[ and used here to compute the spectrum numerically. The quantities allowing us to study 



spectral correlations for periodic quantum systems are introduced in section ffl. Section ^ 
is devoted to the semiclassical theory of the form factor for small times (r ~ '^/N), while 
in section we discuss the form factor for intermediate times (1/A^ ~ r < 1). Finally, in 
section ^, we summarize our results. 

II. QUANTUM MECHANICS 

In the present section we describe the quantization of the system shown in Fig. |l](a). It is 
composed of A^ identical unit cells of length a which form a quasi one- dimensional chain of 
total length L = Na. Each unit cell represents a waveguide of unit width with two half-discs 
as shown in Fig. |^. The radii are denoted by R- and i?+, respectively. The wave function 
■ip{x,y) satisfies the Schrodinger equation 

- 7^— ^^(a;, y) = Eip{x, y) (1) 

with Dirichlet boundary conditions on the walls of the waveguide and periodic boundary 
conditions in x-direction 

ijj{x + L,y) =ilj{x,y). (2) 

Moreover, as shown in Fig. |I](a), the system is invariant under discrete lattice translations 
X — > X + a. In conjunction with (0), the translational symmetry implies the existence of 
discretized energy bands Ea{q). The energies are labeled by a discrete set of Bloch numbers 

g=^, ^ = 0,...,iV-l (3) 

in the first Brillouin zone, corresponding to the N irreducible representations of the group of 
lattice translations. A typical band spectrum Ea{q) is shown in Fig. ^. The wave functions 
can be classified according to their Bloch numbers and can be written as 

^qix, y) = exp(zgx) Uq{x, y) , (4) 



where Uq{x,y) is periodic in x with period a (Bloch's theorem). Hence, for every Bloch 
number q one has the quantization condition 

- ^A ^g{x, y) = E{q) '^jgix, y) , 

il)q{x + a,y) = exp(iga) ^^(x, y) . (5) 

Obviously, the spectrum Ea{q) has the following symmetries 

E^{q + 271/ a) = E^{q), (6) 

E^i-q) = E^iq) . (7) 

In what follows we discuss an implementation of the quantization condition (^. Due to the 
translational invariance, it is sufficient to quantize a single unit cell with periodic boundary 
conditions including an additional phase exp(iga), as shown in Fig. ^ In the straight-channel 
sections of the unit cell, the wave function can be decomposed into normal modes 



oo 



lM fJ+)^+iskiix-x,) , J-)^ 



j=i Jkj 



where 



(j)j{y) = V2 sm{J7Ty), (9) 



and the index s = =F designates the l.h.s. and r.h.s. of the unit cell, respectively. The index 
n labels the normal modes, with wave numbers kn defined according to 



, ^ I y k"^ - (j'tt)^ for open modes (j < A(^)) , . 

iJ{JTTy — k? for closed modes (j > A(£')) , 



where Tik = y2rnE. The number of open modes A(E) is given by the integer part of k/n. For 
an open mode, a^., represents the amplitude of the incoming and a^,, that of the outgoing 
partial wave. As the semicircular obstacles are approached, the partial waves for j > A{E) 
decrease (increase) exponentially. The coefficients of the partial waves to the left and to the 
right of the obstacles can be related to each other using Eq. (|). One has a^-_,_ = e^^°'af/l^ 
which can be written as a matrix equation 



r(^)(r+)j (11) 



%A-^^V^] 



with 

^,.,,v(?) = 5.•,/'^s,-s'e*'''^ (12) 

We will now express the condition (|TT|) in terms of the generalized scattering matrix S{Ei) 



of the obstacle |11] which provides a convenient starting point for both the numerical imple- 



mentation and the semiclassical theory. For this purpose, we consider the Green function 



in the unit cell which satisfies [A + k"^] G{r, r') = —S{r — r'), with Dirichlet boundary con- 
ditions on the walls of the waveguide and outgoing boundary conditions at the borders of 
the unit cell. In the straight-channel sections, G{r, r') can be written in the normal-mode 
decomposition 

3,3 V •? V 3 



+Sjs,rs'e''^^^''-''''^^'''^o'^'''-''''^) . (13) 

The matrix S depends on the particular geometry and remains to be determined. Using (^, 
(p!3|) and Green's theorem for the unit cell it is straightforward to show that S transforms 
the amplitudes of the incoming partial waves into those of the outgoing ones, 

(?)-(:3' 

for any solution of (|l]). Hence, it can be regarded as a generalized scattering matrix for the 
obstacle in the waveguide. The usual unitary S'-matrix is obtained by restricting n and n' 
to open modes. We note that S may be written in block form 

s=(; ,;:), (15) 

in terms of the refiection and transmission matrices r, r' and t, t', respectively. The phase 



shifts are defined as in Ref. |T2[ such that Sjsj's' = Sj,j' Ss,~s' e*'^-'" for an empty unit cell of 



length a. Using ( pTj) to eliminate the coefficients of the outgoing waves in (JT^) we find 

[l-T{q)S{E)] (^r-))=0' (16) 

which can be satisfied by some combination of incoming waves provided that 

det[l-T{q)S{E)] = 0. (17) 

This is the secular equation for the spectrum within the scattering approach to quantization 



rT,TB|. Together with the method described in appendix ^ for the computation of the 
matrix S, Eq. (|l^ can be used to obtain the energy spectrum of the periodic chain of 
scatterers numerically. For this purpose one restricts S* to a finite matrix. In order to obtain 
well- converged eigenvalues, at least four closed modes have to be taken into account in the 
present case. 

III. SPECTRAL CORRELATIONS IN PERIODIC SYSTEMS 

In the present section we discuss the spectral properties of quasi one-dimensional periodic 
systems. A typical spectrum is shown in Fig. ^ (for A^ = 32). The spectral correlations are 
expressed in terms of the density of states 



d{E,q)=Y^5{E-E^{q)). 



(18) 



It is convenient to decompose the density of states into a mean and an oscillatory part, 
d = (d) + d. Note that (d) is the mean density of states for a given Bloch number. In 
other words, it is the mean density of states associated with one unit cell. In order to 
emphasize this, the mean density will henceforth be denoted as (di). As usual, we unfold 
the spectrum [Q in order to obtain a unit total mean density. This is achieved by 

x^iq)=N{d^)E^{q), (19) 

where the energy dependence of the mean density of states has been neglected. Note that 



the unfolding of the spectrum is done with respect to the area of the entire chain [0 We 
define the unfolded density of states according to 



d{x,q) =^5{x-Xa{q)). 

a 

The spectral form factor K^i^q, q'; r) is then defined as 

/ r'^ ~ ~ 

KN{q, q'; r) = { d^ e'^"?^ d{x + ^/2, q) d{x - e/2, q') 



(20) 



(21) 



Here, r is a rescaled time, related to the physical time t by t = 2TTh{di)NT. For a finite 
spectral window the form factor can be expressed as follows [|15 



KN{q,q';T) 



N 
Ax 



X 



E 

a 

E' 



^J2-Ki(x-Xa{q))T 



f{x-Xa{q)) - fir) 



-2m{x-x^,{q'))T f/„_ 



f{x-Xa'{q')) - f{r) 



(22) 



Here, f{x) is a spectral window function defined according to 

fix)- 



1 for |x| < Ax/2 
otherwise . 



(23) 



/(r) denotes the Fourier transform of f{x). Note that f{x) is normalized such that its 
integral equals the number Ax of states contained in the spectral window. 



dx /(x) 



Ax. 



The form factor K^iq-, q''-, t) is normalized such that K^iq, g; r) — * 1 for large r. In fact, for 
large r, only the diagonal contributions to the double sum in Eq. (p2D have to be retained 
and one has (for q and q' not related by the band symmetries discussed in section ^ 

KN{q, q'; t) ^ Sq^q'—J2 f^{x - Xa{q)) 



1 C^ 



5q,q' 



(24) 



The fact that the form factor tends to a constant for q = q' and large r is a consequence 
of the discrete nature of the spectrum. The corresponding quantum dynamics is governed 
by quasi-periodic motion on the time scale of the inverse intra-band spacing. According to 
Eq. (|2^) , correlations between different Bloch numbers decrease rapidly for large r, provided 
the Bloch numbers are not related by symmetries. An example is given in Fig. ^, where 
Ki\f{q, q'] t) is shown for q ^ q' . As expected, the form factor tends to zero for large r. On 
the other hand we observe significant cross correlations between different Bloch numbers at 
intermediate values of r, of the order of r ~ 1/A^. This time scale corresponds to the mean 
inter-band spacing and the correlations are associated with the existence of energy bands. 
Before discussing these cross correlations in detail, we consider the spectral correlations at 
fixed Bloch number q. 

A. Spectral correlations at fixed Bloch number q 

In the present subsection we discuss spectral correlations for fixed Bloch number q. Consider 
for instance the spectrum shown in Fig. 0. Since the classical dynamics of the unit cell is 
chaotic, the spectral correlations at fixed q are expected to be described by random-matrix 
theory. The quantization condition (^ of section || can be rewritten as follows 



2 



h 

2m 



Uq{x + a,y) =Uq{x,y). 



Obviously, the Bloch number q can be regarded as reflecting an Aharonov-Bohm flux 0, 
with q = 27r/a 0/0o and 0o = 2iTh/e. For g = 0, the Hamiltonian H{q) is invariant 
with respect to time reversal and the spectral properties are described by the ensemble 
of Gaussian orthogonal random matrices (GOE). In the case q = vr/a, H{q) is invariant 
under the anti-unitary transformation 6 = exp{—2iqx)J, where J is the operator of complex 
conjugation. For q = vr/a, the appropriate ensemble is hence also the orthogonal one (GOE). 
For q = 7r/2a, on the other hand, the spectral properties are described by the ensemble of 
Gaussian unitary random matrices (GUE). As q changes from zero to it /2a, one expects a 
crossover from orthogonal to unitary statistics. 

Fig. ^ shows the spectral form factor i^7v(g, g;r) as a function of r for two different 
values, g = and q = 7r/2a, of the Bloch number. The universal results for the orthogonal 
and unitary ensembles [|I^ are displayed with solid and dashed lines, respectively. As ex- 
pected, we observe GOE behaviour for q = and GUE behaviour for q = 7r/2a. Similarly, 
the histogram -P(A) of nearest-neighbour spacings exhibits a crossover from GOE to GUE 
statistics. This is shown in Fig. |^. This crossover occurs on an ^-dependent scale. In the 
semiclassical limit h —>■ 0, the transition is instantaneous. 



B. Spectral cross correlations 

In the present subsection we discuss correlations between different Bloch numbers. These 
correlations decay rapidly for large r, as shown for instance in Fig. ^. However, it will 
turn out in the following that spectral cross correlations at intermediate values of r are 



of particular importance in extended periodic quantum systems, since they are directly 
related to the existence of discretized energy bands. In order to analyse correlations between 
different values of g, we define 

S^{n; r) = ^Y. e^"('^-^>ir^(g, q'; r) . (25) 

q,q' 

Note that S'Ar(0;r) is just the conventional spectral form factor Kn{t) for the composite 
system referred to in the introduction. Rewriting 5'iv(n;r) in terms of a Bloch-number 
specific density of states 

d{x,n)=J2e"''"'d{x,q), (26) 

<? 

we obtain 

SN{n;T) = I r di e~^^'^^ d{x + i/2,n)d{x - il2,ny) . (27) 

In the following, we discuss the long-time asymptote of Sn^u-jt). For a finite spectral 
window, SnIu; t) can be expressed in a form similar to Eq. (^) and for r ^ 1/Ax one has 

2 



Mn;r)^i- 



(28) 



q a 

Taking into account the symmetry of the spectrum with respect to g ^ 27r/a — q and in the 
absence of further symmetries, the phases of the off-diagonal terms in the double sum are 
random and cancel. Using q = 2ttv/ L, one has 



1 

iV 



N/2-1 r, 

4-^ r, /ZTnyn 



1 



2 + 4 gcos^ (_)]-/_ J. m (29) 



assuming A^ to be an even integer. Using Jdx f'^{x) = Ax (compare Eq. (P5|)), we obtain 

for n = 0, N/2 , 




(30) 

For large A^, Sisfin; r) — > 2 for n = 0, N/2. This is a consequence of the band symmetry (^ 
discussed in section y. 

IV. SEMICLASSICAL ANALYSIS 

In the following we obtain a semiclassical approximation for Sj\f{n; r). This is done in two 
steps. First, we derive a semiclassical approximation for the density of states for a periodic 
quantum system. It will be seen below that in the present case it is particularly convenient 
to start from Eq. (|l^. As discussed in Ref. 0, the resulting semiclassical expression features 
the characters x(^5 q) = exp(mga) of the group of lattice translations. Second, we introduce 
the diagonal approximation for the form factor. The diagonal approximation is adequate 
for times smaller than the Heisenberg time of the unit cell, that is for t r^ tu = 27rh{di) 
or equivalently for r ~ 1/A^. A theory covering the long-time behaviour is worked out in 
section 0. 



A. Semiclassical approximation for the density of states 

We now proceed to derive a semiclassical approximation for the energy spectrum. The start- 
ing point is Eq. ([T7|), in which S{E) is restricted to the unitary A x A-matrix corresponding 
to the open modes. In this case the spectral density following from (|l^ may be written as 



11 



diE,q) = Y.6iE-E^iq)) 



a=l 



A^^ + ^^Im E — Tr [T{q)SiE)Y 



(31) 



where Q{E) = —i log det S{E) denotes the total phase of the S'-matrix. Using a theorem by 
M. G. Krein [|17|, it is possible to obtain a semiclassical approximation for B which contains 
a smooth part given by the generalized Weyl law W^ for the unit cell and an oscillating part 
which represents the standard Gutzwiller contributions from all the periodic orbits which 
are trapped inside the unit cell |]TU[. In order to obtain a semiclassical approximation for 
Tr (rS*)™, we first use (0) at x = Xs, x' = Xg' and find 



i^S)js,j' 



2s e 



isqa 






(32) 



The completeness of the normal modes Y.'jLi 4>j{y)4>j{y') = ^{u ~ u') for < y, y' < 1 allows 
to express the powers of TS in coordinate space as multiple integrals over the Green function 



± 
Tr (r^)"* = 2 V^'^ E fdyi...dy„ 



X 
X 



si-g—G{xi,yi;Xs2,y2) 
d 



X\=X- 



dXr. 



' ^yXrifYi^ i/mi '''su i/lj 



where the winding number n is defined as 

— n = si + V Sjn- 

The standard semiclassical expression for the Green function 

-1/2 



G{x,y]x\y') = -^ 



27ik 



dx I 



dip' 



exp ikL 



71 



(33) 



(34) 



(35) 



can now be inserted into (p3D. The sum in (|35D runs over all classical trajectories a leading 
from {x',y') to {x,y). L^ is the length of the path a and /i^ contains the Maslov index of 
the trajectory plus twice the number of refiections from the walls of the waveguide. Due 
to the special geometry the Maslov index happens to be zero for all trajectories, ip and x± 
denote the direction of the trajectory and the coordinate normal to it, respectively. 

The integrals in ( P^D are evaluated in the saddle-point approximation. The condition of 
stationary phase is satisfied for all periodic orbits p of the unit cell with periodic boundary 
conditions, i.e., whenever a periodic orbit leaves the unit cell to the left or to the right, it is 



9 



reinjected from the opposite side. Tr (TS)"^ contains all those orbits for which this happens 
exactly Up = m times, irrespective of which of the two boundaries is hit. In contrast, Up 
defined in ( p^ ) counts the number of reinjections including a sign according to the direction, 
such that Up finally represents the number of windings of the periodic orbit around the unit 
cell (see Fig. |^. Obviously the orbits contributing to the total phase have a winding number 
Hp = and for the other orbits we have |?T-p| < rrip. We obtain 

d{E, q) = — — Y^ Wp^rTp exp (^^Sp + irUpqa - i-rfipj , (36) 

where r = — cx), . . . , cxo (excluding r = 0), p is summed over the primitive periodic orbits, 
u!p,r = I det(Mp — 1)|~^/^, M = d{xj_, (f)/d{x'j_, ip') is the stability matrix of the periodic orbit 
p, and Sp = hkLp is the corresponding action. The discrete translation symmetry results in 



an additional phase which depends on the winding number of the orbit. Eq. (|36D represents 
a special case of the general result obtained in ^. According to Ref. |]^, the additional 
phase factors are related to the characters x(^5 q) of the irreducible representations of the 
group of lattice translation. 



B. Diagonal approximation for S'Ar(n, r) 

In the following we derive the semiclassical approximation for S]\f{n;T) in the diagonal 
approximation ^^. Using Eqs. (|19|) and (^), we have from ( p6| ) 

d(x,n) ''"^ " ~ 



2^^?^ 5 ^'rUrl Wp,rTp exp ^-rSp + r-r^^p 



Here 6^^ denotes a Kronecker delta with the arguments taken modulo A^. Expanding the 
actions according to Sp{E ± e/2) ~ Sp{E) ± Tpe/2, we have 

S^{n- t) = ( r d^ e~^"^«" d{x + ^/2, n)d{x - ^/2, n)) 



{2'Kh{d 






'^p.r'^p' y -^p ^p' 



f % % TV 

X exp [TirSp - r'Sp') + ^{rTp + r'Tp,)e + i-{riip - r' iip, 



where ^ = N{di)e. Provided the actions Sp are uncorrelated, the off-diagonal terms in 
the double sum over primitive periodic orbits vanish upon averaging. This is certainly the 
case for periodic orbits with periods less than some fraction of the Heisenberg time of the 
unit cell, Tp r^ tn = 27r^(o?i), which restricts the applicability of the present approximation 
to times smaller than the Heisenberg time of the unit cell, or equivalently, to r ~ 1/A^. 
Furthermore we note that pairs of time-reversed orbits have in general opposite winding 
numbers. However, periodic orbits with np = or A^/2 give rise to an additional degeneracy, 
since in these cases —Up = np (modA^). We define the coherence factor 



10 



2 forn = 0,iV/2, 

1 1 otherwise , 

and obtain 

S^in; t) ~ 7„riV2 ^ ^f^)^^ |u;p_,|2 ^^ ^ ^^^^ _ 27r;i(di)iVr) . (38) 

p,r 

For n = 0, the form factor S'Ar(0; r) is just the conventional spectral form factor Kn{t) for 
the composite system, i.e., for a billiard chain with periodic boundary conditions (0). The 
primitive periodic orbits of the entire chain have winding numbers n^ which are multiples 



of A^, as opposed to the periodic orbits of the unit cell. The Kronecker delta in Eq. (38 
ensures that for n = 0, only periodic orbits of the entire chain contribute. 

The sum over periodic orbits in Eq. (|3^) can be evaluated as follows. According to 
appendix C, it can be replaced by the coarse-grained classical propagator P(n; t) defined in 
Eqs. ( pi|) and ( piO|) . Furthermore, the classical dynamics in the periodic chain is diffusive 
in the x-direction (provided direct trajectories are eliminated). This is due to the chaotic 
scattering in the unit cells. Under these circumstances, P(?t,; t) is given by the expression 
( pSf ) derived in appendix C In order to show that the classical dynamics in the billiard 
chain is adequatly described by Eq. ([B5|) , we have performed numerical simulations of the 
classical dynamics in periodic chains. The chains consist of iV = 8 or A^ = 128 copies of the 
unit cell shown in Fig. ^ The radii of the half discs have been chosen large enough so as to 
eliminate direct bouncing-ball trajectories {R- = 0.8 and -R+ = 0.4 in units of the channel 
width). The lengths of the straight-channel sections are ~ 0.1 in the same units. The results 
of the numerical simulations are shown in Fig. ^ (solid dots), as a function of the trajectory 
length / = hkt/m. For a given trajectory length /, we have started 10^ trajectories and 
calculated the average probability to move n unit cells away from the starting point. For 
convenience, the results are given as a function of lengths / rather than times t. Also shown 
is a fit to Eq. (|B5|). The fitting parameter is the diffusion constant D. In the case shown, 
we obtain for the scaled diffusion constant Dm/hka^ ^ 0.082. We observe good agreement 
between the numerical simulations and the results of Eq. (P5|). 

Replacing the sum over periodic orbits in Eq. (|38D by the appropiate expression for the 
coarse-grained classical propagator, we obtain (neglecting repetitions of periodic orbits, since 
primitive periodic orbits proliferate exponentially) 

SN{n; t) ~ -fnTN'^P{n; 27rhN{di)T) 

where we have introduced the dimensionless parameter 

2TT'^h{di)D 



9i 



a? 



The parameter gi is analogous to the dimensionless conductance of disordered conduc- 
tors [|l^. Note that it is defined with respect to the unit cell. However, since we do not 



assume static disorder within the unit cell, the analogy is purely formal. 
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Eq. (|39|) is the central result of this section. It relates the quantum spectral form factor 
S]\f{n;T) to the classical diffusion propagator P{n;t) of winding numbers n. According 
to Eq. (|39|), the behaviour of Siy{n;T) depends crucially on the values of A^ and gi. We 
distinguish two cases, gi/N"^ > 1 and gi/N'^ < 1. The spectral properties are qualitatively 
different in these two regimes. For gi/N'^ > 1, the classical dynamics becomes ergodic 
before the energy-time uncertainty relation allows to resolve the inter-band spacing. This 
corresponds to a time scale of the order of the Heisenberg time of the unit cell, tn = 2TTh{di), 
or equivalently r ~ 1/A^. The sampling of the energy bands by discrete Bloch numbers q is 
then too coarse to reveal that continuous bands are underlying the discrete levels, and the 
full spectrum appears as a superposition of A^ independent spectra. 

For gi/N"^ < 1, on the other hand, we have giT < N throughout the range of validity 
of Eq. (|39|) (0 < r ~ 1/A^) and the sum in (|39|) can be restricted to /i = 0. Classically 
this means that the diffusion cloud generated by the propagator of winding numbers is still 
well localized at times corresponding to the Heisenberg time of the unit cell (r ~ 1/A^). 
Semiclassically, the coarse-grained diffusion propagator encodes the quantum-mechanical 
spectral correlations due to the existence of energy bands. At times of the order of r ~ 1/A^, 
where individual bands begin to be resolved, the spectrum appears to be A^-fold degenerate. 
As anticipated in the introduction, the clustering of the levels into bands is reflected in 
the spectral form factor: at r ~ 1/A^, the form factor exhibits a peak of the order of 
A^, Sn{0]t) ~ A^. Hence, for gi/N^ < 1, the semiclassical expression Eq. (pQ]) describes 
the effect of the quantum levels clustering into bands and relates this to the space-time 
dependence of the classical diffusion propagator. 

In the remainder of this section, we discuss the signatures of level clustering in numerical 
band spectra and compare with the predictions of Eq. (^). We have calculated the quantum- 
mechanical form factors Si\f{n; r) for chains with N = 8 and A^ = 128 unit cells, respectively, 
and with the same geometrical parameters as used in the classical simulations described 
above. Since S^in; r) is not self-averaging with respect to energy, we have taken the mean 
of seven samples with slightly different geometries (the relative deviations are of the order 
of a few de Broglie wavelengths). In addition, we have performed an energy average, around 
a mean energy of 2mE/h'^ ~ 8.9 x 10^. For unit channel width, this corresponds to 30 open 
channels. Fig. ^(a) shows the quantum data for A^ = 128 and winding numbers n = 0, 2, 4, 8 
(solid lines). Fig. ^(b) shows the quantum data for A^ = 8 and winding numbers n = 
0,1,2,3 (solid lines). In both cases, the semiclassical approximation according to Eq. ( PU[ ) 
is also shown (dashed lines). We observe reasonable agreement with the semiclassical result 
for times r ~ 1/A^. For larger times, the diagonal approximation leading to Eq. (|39D is 
inappropriate as remarked above. In the present section, we restrict our discussion to times 
r ~ l/N. 

For A^ = 128 as well as for A^ = 8, the quantum data for zero winding number [n = 0) 
are enhanced by a factor of two, in keeping with the prediction of Eqs. ( p^ ) and (p9|). For 
n = and for very small times r <^ 1/A^, we observe a square-root dependence of the form 
factor, 5*^^(0; r) ~ i/r. This law persists up to r ~ 1/A^ and can be understood in terms of 
the behaviour of P(0;t). For small times, P(0; t) ~ l/Vt- In Eq. (^) this corresponds to 
keeping only the /i = term which leads to Sn{0;t) ~ y/r. For very short chains, on the 
other hand, corrections due to the finite length modify this behaviour when the n ^ terms 
contribute significantly. 
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In both cases, the conductance is gi ~ 33. In the first case (A^ = 128) we have gi/N"^ < 1 
and the structure of Sn{0; r) at r ~ 1/A^ (see Fig. ^(a)) is the signature of levels clustering 
into bands, as discussed above. As gi/N"^ grows larger, on the other hand, the effect of 
clustering of levels disappears (see Fig. ^(b), where N = 8 and hence gi/N"^ — 1). The 
sampling of the bands by discrete Bloch numbers is then too coarse to reveal the underlying 
continuous bands. 

In summary we emphasize that for very small times the quantum dynamics of extended 
periodic quantum systems is governed by diffusion. For small gi/N"^, the spectral form factor 
exhibits a peak ~ A^ associated with the clustering of levels into bands. 

To conclude we mention a sum rule for the spectral form factors. According to ap- 
pendix C, the diffusion propagator P{n; t) is normalized to J2n=o P{f^] t) — 1- Using 

1 N-l N-l 

^ E 7n-'^7v(n;r) ^ rA^ ^ P{n-2Tih{d^)NT), 

-'^ n=0 n=0 

this implies a sum rule for the spectral correlators SN{n; r) for r < 1/A^. Together with (|30D 
we have for large A^ 

iv5/" M-;r)-|^ ^^^^^^/_ (40) 

This sum rule is shown in Fig. |1^. Also shown is the expectation according to Eq. (^). 

V. BALLISTIC SPREADING AND THE SCALING OF THE FORM FACTOR 



In sections |T| and |I^, the behaviour of the spectral form factor Sn (n; r) has been determined 
for very large times and for small times. In the first case, for times of the order of the inverse 
intra-band spacing, we have 5'Ar(0;r) ~ 7^, corresponding to quasi-periodic motion. In the 
second case, for r ~ 1/A^, the quantum dynamics is diffusive. In this regime, and for 
gi/N"^ < 1, spectral correlations due to the clustering of levels lead to a pronounced peak in 
Sn{0] r), at times of the order of r ~ 1/A^. How are these two regimes to be connected? In 
the present section we develop a theory for the spectral form factor at intermediate times. 
The derivation of Eq. (^) is based on the diagonal approximation which is known to fail 
for times r > 1/A^. In the limit of large A^, on the other hand, the sum over q in Eq. (p8[) 
can be replaced by an integral 

L 



SN{n;T) ~ — 



— fdq E e*"''"-2™.{g)ry^^ _ ^^(^^^ 



(41) 



According to ( |19|) we have 27rxQ,(g)r ~ Nt and hence the phase of the integrand is rapidly 
oscillating for large r. In this case, the integral can be evaluated in the saddle-point ap- 
proximation. The condition of stationary phase reads 2nx'^{qj) = na/r or, in terms of the 
unsealed energies Ea{q), 

""" na 



Q pi 

2n{d^)N ' 



dq 



(42) 
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We emphasize that the saddle points qj depend on the parameters r, n and A^ only through 
the combination n/Nr. The resulting expression for 5'Ar(n;r) is 



SN{n;r)^^ 



2^ "ir v^^fej^ 



The sum over a and Qj involves many terms with uncorrelated phases. Assuming some 
smoothing over r, we can justify replacing it by 

Mn;r) = ^F(|-), (43) 

where 

\NtJ An^Ax^^ {d,)E'^iq,) ' 

In writing (^3|), we have used the fact that the saddle points are functions of n/Nr. For 
n = 0, in particular, Eq. (|^) implies Sn{0; t) ~ 1/r. According to Eq. (^21), this behaviour 



corresponds to a ballistic quantum dynamics, J {'n?) ~ r. This is in contrast to the diffusive 
spreading derived for short times, {p?) ~ r for small r. 

Clearly, the ballistic spreading is limited by the finite system size. We denote the time it 
takes for the distribution in n to spread over the entire system by r*. This is a time of the 
order of the inverse intra-band spacing. According to the discussion in section |Tl|, r* is just 
the time where the form factor S]^{n] r) saturates at its asymptotic value. This saturation. 



discussed above in section [1IIB| , is not described by Eq. (^31). 

For n = 0, Eq. ( ^Hj ) takes a particularly simple form. Matching (^) at r = 1/A^ with 
the peak derived in the diagonal approximation, we find 

F(0) = l/y2^. (44) 

Furthermore, since S'7v(0;r*) ^ 2, we have r* ~ F(0)/2. The behaviour of S'Ar(0;r) for 
gi/N"^ < 1 is now completely determined on all time scales. We have 

r F(0) iVVwV for < r < l/N 
Sn{^] r) ~ 2 J F(0)r-i for l/N ^t< F(0)/2 (45) 

[l forr>F(0)/2. 

This prediction is tested numerically in Fig. |TT[ 5'Ar(0; r) is shown for various values of N 
with otherwise unchanged parameters in (a) and for one particular value of A^ on a doubly 
logarithmic scale in (b). The dashed lines show the prediction of (^5]). Fig. 0(a) shows that 
the A^-dependence of the quantum data are adequately described by Eq. (^). In Fig. pl](b), 
one can distinguish several time regimes, as expected from Eq. (^): First, for very short 
times the form factor increases according to ~ ^/r, as predicted by Eq. (^), which reflects 
the underlying classical diffusion. This behaviour continues up to r ~ 1/A^. Second, the 
peak in SjqiS^] r) at r ~ 1/A^ is a signature of the levels clustering into bands. Accordingly, 
the height of the peak scales as ~ A^. Third, for larger times, 5'Ar(0; r) decays as ~ 1/r. This 



is consistent with a ballistically spreading quantum distribution, J {th?) ~ t. In this regime 
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the curves corresponding to different values of A^ coincide (see Fig. 0(a)) as predicted by 
Eq. (^5|). Fourth, the form factor saturates at r*. 

As A^ — > oo, the spectral form factor approaches a universal function. The diffusive do- 
main < r < 1/A^ shrinks and for t < t* the form factor is dominated by the 1/r behaviour 
which reflects the formation of continuous bands. This behaviour was not discussed before, 
and is one of the most important results of the present work. 

We should finally note that periodic systems in d dimensions will exhibit a t~'^ behaviour, 
as can be easily seen from the derivation above. 

VI. CONCLUSIONS 

In the present article, we have investigated the implications of a discrete translation sym- 
metry for the spectral correlations of quasi one- dimensional quantum systems with chaotic 
classical dynamics. The decisive new feature here, as compared to disordered systems, is the 
presence of discretized energy bands. We have introduced a generalized spectral form factor 
involving spectral correlations both for identical as well as for non-identical Bloch numbers. 
We have shown in detail how the time-dependence of this form factor reflects the quantum 
dynamics in the periodic system: 

(i) For small times, the form factor was shown to be semiclassically equivalent to the 
coarse-grained classical propagator. In the present case, chaotic scattering implies classical 
diffusion in the extended system. This enabled us to derive an explicit expression for the 
form factor, depending only on the diffusion constant and the system size. 

(ii) At r ~ 1/A^, the spectral form factor exhibits a peak due to level clustering provided 
the discretization by the Bloch numbers is sufficiently fine in order to reveal the underlying 
smooth bands. The height of the peak corresponds to the number of strongly correlated 
levels and hence it is proportional to A^. 

(iii) For 1/A^ ~ r < r*, we have derived a universal scaling law for the form factor. This 
regime corresponds to ballistic quantum dynamics and is restricted by the time r* needed 
to ballistically cover the entire system. 

(iv) For T > T*, the form factor attains its asymptotic constant value. This regime 
corresponds to quasi-periodic motion. 

We emphasize once again that our results — tested numerically using the example of 
a quasi one- dimensional periodic chain of chaotic billiards — are more generally valid (see 
also |T^) and in fact universal. Finally we note that the transition from a periodic chain 
to a quasi one-dimensional disordered system is also governed by a universal law. A study 
concerned with the effect of introducing disorder in the present framework is under way. 
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APPENDIX A: THE 5-MATRIX OF THE UNIT CELL 

In this appendix we briefly describe the numerical computation of the generalized S'-matrix 
of the unit cell needed according to ([T7|) for the quantization of the periodic billiard chain. 
As shown in Fig. ^, the unit cell consists of two half discs with radii Ri and i?2- For 
the method we employ here, it is essential that the projections of the half discs onto the 
X-axis do not overlap (see Fig. ^. This allows us to express the total S'-matrix in terms 
of the reflection and transmission matrices of a single half disc in a waveguide using the 
concatenation formulae 



t = t2{l + r[[l - r2r[]-^r2)h 
r = ri + t'^[l - r2r'J"V2ti 



t' = t\[l - r2r[]-% 



(Al) 



2 



r' = r2 + t2r[[l — r2r[] ^t'2 



2 



Here ri, r[ and ti, t'^ describe the scattering off the left half disc, while r2, r'2 and ^2, t'2 describe 
the scattering off the right half disc. The definitions of these matrices are analogous to the 



definition of r, r' and t,t' in (JT^) and (|T^). 

The computation of the S'-matrix of a single half disc inside the waveguide remains to 
be explained. We assume that its center is located at the origin. Other reference points 
as well as the reflection from y = 1/2 (which is needed to obtain the right obstacle in 
Fig. H) can be obtained in terms of unitary transformations. Moreover we use the reflection 
symmetry of the half disc in order to decompose the S'-matrix into a symmetric part S"*" 
and an antisymmetric part S~, which are related to the reflection and transmission matrices 
according to S^ = t ± r. The calculation of S~ including evanescent modes has been 
described in detail in |TB[ and S~^ can be found by straightforward analogy. The result can 



be summarized in the matrix equation 

S± = l + C±T(P±)-^(C±)*^ (A2) 

Here T denotes the transition matrix for the scattering from a circle in the free plane 
which is diagonal in angular momentum representation and has the elements Ti{kR) = 
-Ji{kR)/H^{kR). The matrix 

P± = 6u' - [Fi-v{2k) T Fi+i\2k)]Tv{kR) (A3) 

with 

F,(a;) = 2 5:cos(^-jif+(na;) (A4) 

describes the coupling of angular momenta due to the presence of the waveguide. The 
angular momentum runs over all positive even numbers for S~ and over all positive odd 
numbers for S^ . For the numerical computation of the only conditionally convergent series 
(|A^), a method has been developed in |^. Finally, the matrix 




Cti = isi-(Ru^R^.-i) (A5) 
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with 



it,- 



— 2Sgn(/) kj + jtt' 
k 



(A6) 



contains the expansion coefficients of the normal mode n in the free-channel region in the 
angular-momentum basis which is appropriate close to the semicircular obstacles. The 
method described in this appendix has been successfully applied to compute S~ up to 
values of fc ~ 1000 |jl9|, which is far beyond what is needed within the present application. 
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APPENDIX B: CLASSICAL MECHANICS OF PERIODIC CHAINS 



In this appendix, we discuss the classical mechanics of a periodic chain of chaotic scattering 
elements with periodic boundary conditions. We assume that the chain extends along the 
x-direction and consists of A^ identical copies of a unit cell (of length a). Let z = {x,y,ip) 
denote a phase space vector on the energy shell {y and ip are the transversal coordinate and 
the direction of the momentum, respectively). The dynamics on the energy shell is governed 
by the propagator 



P(z,z';t) =6(Z{z',t)-z) 



(Bl) 



where Z(z',t) denotes the trajectory started at t = from z'. Since the unit cell is chaotic 
(neglecting the possible influence of an infinite horizon), both y and ip are effectively decor- 
related on time scales larger than the ergodic time of the unit cell. Consequently, we obtain 
diffusive spreading in the x-direction 



The classical propagator solves the diffusion equation 

P{x,x';t) = 



dt 2 9x2 



(B2) 



with the boundary condition P(x; 0) = 5^^\x), where 5^^\x) is the L-periodic delta function. 
The propagator is normalized according to /q dx P{x]t) = 1. The solution of (|B^ ) can be 
easily found by representing P{x, t) in terms of eigenf unctions of the stationary problem. 
We obtain 



P{x;t) 



L 



J2 exp 



2T(iu- 



X 



D 



27rz/\2 t 



(B3) 



L \ L J 2 

It is appropriate to introduce an additional coarse graining of the coordinate x by defining 

[x/a] (B4) 



n 



where [. . .] denotes the integer part. As will be seen below, n can be interpreted as the 
winding number of a trajectory around the unit cell (with periodic boundary conditions). 
^From Eq. (|B3| ) we have 



/ N 1 ^ ( n D 
Pin; t) c::^ — > exp 2^^ 

which after Poisson summation can be written as 



27rz/ 



n2 



N 



t 



P{n;t) 



^2tx Dt 



Y. 6^P 



2Dt 



[n - fiNf 



(B5) 



(B6) 



Due to the finite system size, the propagator P{n; t) exhibits a crossover from diffusion to 
equidistribution on the time scale of the Thouless time to = L'^/'n'D. Indeed for small times 
t — > 0, only the /x = term contributes in (|B^ ) and we have 



n^a^ 



Pin] t) ~ , exp . (B7) 

^ ' V2^Fm ^ \ 2DtJ ^ ' 

For large times (t ^ Id) the dominant term in ( P3|) is i^ = and P{n\ t) saturates at 

P(n;t)^i^. (B8) 



This behaviour is summarized in Figs. |T2| and |T3[ Fig. |I2| shows the diffusion of winding 
numbers. The solid dots are a numerical simulation for an infinite periodic chain. The solid 
line is a fit to the diffusion law (x^) = Dt, or equivalently (n^) = Dm/hko? I. We obtain 



Dm/Uka? ~ 0.082. Fig. |T3| shows P(n; t) as a function of n for three cases: (a) for t ^ ^d, 



(b) for t c:^tD and (c) for t > t-Q. 

So far we have derived a diffusion equation for the coarse-grained propagator P{n\ t) by 
neglecting all details of the chaotic dynamics within the unit cell. In the following we derive 
an alternative expression for P{n\ t) in terms of all periodic orbits of the unit cell, which 
will later allow to make contact with the semiclassical spectral density. For this purpose we 
employ (pl|) and ( p^) to write the coarse-grained propagator as 

P{n- 1) = — ^ / dz dz'P(z + na, z'; t) , (B9) 

where a = {a,y = 0,ip = 0) and fli{E) denotes the volume of the energy shell for motion 
within the unit cell (as opposed to the volume of the energy shell for the entire chain). 
Since the dynamics within the unit cell is chaotic, the transition probability is effectively 
independent of z', and we use this to replace z' by z: 

^(^' ^) - TTTi^ / dzdz'P(z + na,z;t) 

= f dzP(z + na,z;t). (BIO) 

JUiiE) 

We have now expressed the coarse-grained popagator by the probability to return to the 
same phase-space point up to a translation by n lattice vectors or, when interpreted in the 
unit cell, as the joint probability to return and to have completed n windings around the 
cylinder. The total probability to return on a compact phase space can be expressed as a 
sum over all periodic orbits [l^ and an analogous relation for P{n; t) can be obtained in 
terms of the periodic orbits of the unit cell with winding number n (disregarding bouncing 
ball orbits). We insert (|B1| ) into (|B10|) and use a linear approximation for Z(z,t) in the 



vicinity of each periodic orbit p. x and y are replaced by local coordinates x± and x^ aligned 
with the orbit, where the latter is cyclic with the length Lp of the orbit x\\ +Lp = x\\. Finally 
we have 

P{n; t) ^ E / dx|| dx^dip 6^%^ 5^'^^\Lp + 2k[t - vTp]) 

JdxAt), dxxit) ^ . \ J difit) . dip(t) ^ . \ 



p,r I ^^''v^'-'p 
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where Tp denotes the period of the periodic orbit and Mp is the monodromy matrix. The sum 
over r is a sum over repetitions of primitive periodic orbits p. In this way we have expressed 
the coarse-grained classical propagator of the extended billiard P{n,t) using periodic orbits 
of the unit cell. We make use of this result in the semiclassical analysis of the billiard in 
section ^\. 
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FIGURES 
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(b) 




FIG. 1. Quasi one-dimensional chains of coupled billiards, (a) shows a periodic chain while 
(b) shows a disordered chain. In both cases periodic boundary conditions are imposed at the ends 
of the chains. 
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FIG. 2. Typical band spectrum of a periodic billiard chain with A^ = 32 unit cells. Shown are 
the levels ka{q) as a function of the Bloch number q in the first Brillouin zone, —it /a < q < n/a. 
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FIG. 3. Unit cell of the periodic chain of chaotic billiards. The unit cell has unit width and a 
length a. According to Bloch's theorem, the periodic chain can be quantized using Eq. (p!7|). This 
corresponds to quantizing a single unit cell with periodic boundary conditions and an additional 
phase exp(iqa). 
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FIG. 4. Spectral form factor K]\f(q,q';T) for q ^ q' as a function of tN. For large times 
T > 1/A^, Kiy(q, q'; r) tends to zero since levels of the same band decorrelate for q ^ q' . The arrow 
indicates r ^ 1/iV, a time of the order of the inverse inter-band spacing. The data are taken from 
a chain with N = 128 unit cells. 
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FIG. 5. Form factors K]sf{q,q;T) for different values of the Bloch number g as a function of 
tN. In (a) we show the case q = where the spectrum shows GOE statistics, and in (b) the case 
q = 7r/2a where the spectrum shows GUE statistics. 
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FIG. 6. Histograms -P(A) of the nearest-neighbour spacings A, for fixed Bloch number q. (a) 
shows the case q = where the spectrum shows GOE statistics, while (b) shows the case q = 7r/2a, 
where the spectrum shows GUE statistics. 
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FIG. 7. Topology of periodic orbits in the periodic chain, (a) shows a periodic orbit with 
winding number n = 1 while (b) shows a periodic orbit with winding number n = 0. Note that 
the periodic orbit in (b) is self-retracing. 
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FIG. 8. Diffusion propagator P{n; t) for a periodic chain, (a) shows the distribution of P(n; t) 
for A^ = 128, (b) shows P{n; t) for A^ = 8, as a function of I = fikt/m (the energy dependence of the 
classical mechanics has been scaled out). Note that for large i, P{n;t) approaches 1/A^ (compare 
Eq. (il)). 
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FIG. 9. Shows S]s[{n;T) as extracted from the quantum spectrum for a periodic chain, (a) for 
N = 128 and n = 0, 2, 4, 8 and (b) for TV = 8 and n = 0, 1, 2, 3 (sohd hues). Also shown (dashed 
lines) are the results of the semiclassical approximation Eq. (| 
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FIG. 10. Semiclassical sum rule. The solid line is the sum of the spectral form factors S]\f{n; r), 
while the dashed line is the semiclassical prediction according to Eq. (|40|). 
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FIG. 11. Scaling of S'Ar(n;r). (a) shows the scahng of S'iv(n;T) with N, for A^ = 8, 16,32,64 
and N = 128. (b) shows the various regimes in a doubly logarithmic plot: (i) the diffusive increase 
~ "v/r for small times, (ii) the peak ~ A^ at r = 1/A^ due to the level clustering, (iii) the 1/r decay 
associated with ballistic quantum spreading and finally (iv) the constant asymptote at the time 
scale corresponding to the mean intra-band spacing due to quasi-periodic motion. 
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FIG. 12. Diffusion in the periodic chain. Shown is the the result of a classical simulation: 
The variance {v?) of winding numbers n for trajectory segments of length I increases according to 

(n2) ~ /. 
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FIG. 13. The distribution of P{n; t) for three cases: (a) for t ^ to, (b) for t — to and (c) for 



t > to- The dashed Hnes show the hmiting distributions ( |B7[) and (pq) , for small and large times, 
respectively. 
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